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Molecules at the air-water interface often form inhomogeneous layers in which domains of different 
densities are separated by sharp interfaces. Complex interfacial pattern formation may occur through 
the competition of short- and long-range forces acting within the monolayer. The overdamped 
hydrodynamics of such interfacial motion is treated here in a general manner that accounts for 
dissipation both within the monolayer and in the subfiuid. Previous results on the linear stability 
of interfaces are recovered and extended, and a formulation applicable to the nonlinear regime is 
developed. A simplified dynamical law valid when dissipation in the monolayer itself is negligible is 
also proposed. Throughout the analysis, special attention is paid to the dependence of the dynamical 
behavior on a characteristic length scale set by the ratio of the viscosities in the monolayer and in 
the subphase. 



I. INTRODUCTION AND EXPERIMENTAL 
BACKGROUND 

Molecular layers of surfactants or polymers at the air- 
water interface are often found in inhomogeneous states 
within which appear domains of nearly uniform den- 
sity KB. In many ways, these states resemble con- 
ventional two-phase coexistence regions with sharp in- 
terfaces between different homogeneous phases; because 
of long-range electrostatic interactions within the mono- 
layer, however, domains of a given phase can be stable 
rather than coarsening in time ||. Experimentally ob- 
served domains typically have sizes of 10 — 100 //m. 

A considerable body of experimental [p|-^2| and the- 
oretical ]l3|-|l5|] work has focused on the motion of the 
domain boundaries. These investigations have had two 
complementary motivations. On the one hand, the 
boundary dynamics provide a means of probing physical 
parameters, such as the line tension between phases, that 
are otherwise difficult to measure 0,|l^,|l5|. On the other, 
the electrostatic interactions are caused by the molecules' 
permanent dipoles that are oriented with respect to sur- 
face of the water; their mutual repulsion can result in 
intricate fingering instabilities that have parallels in a 
variety of other pattern- forming systems |3|]. The laws 
of motion of the interfaces between monolayer phases 
have been probed directly in experiments monitoring 
the relaxation of domains to a circular ground state, 
starting either from an elongated "bola" shape @,|l^] or 
from smaller elliptical deformations j^jjjfl- The fastest- 
growing mode at the onset of a fingering instability has 
also been examined P.f)|jTT| . 

The complete interpretation of such experiments re- 
quires an understanding of the hydrodynamics of thin 
layers coupled to a subfiuid. Some progress has been 



made in the analysis of particular cases. Building on 
earlier work on the diffusion of a cylinder imbedded 
in a membrane pi, Stone and McConnell have solved 
the linearized boundary dynamics about a circle p3| , |l4| . 
Schwartz, Knobler, and Bruinsma |l7]] and Stone |Tq ] 
have examined the flow of monolayers through chan- 
nels. The related problem of fingering in quasi- two- 
dimensional domains of ferrofluids has also recently been 
considered ]l9|-p3| , as has the behavior of capillary waves 
in the presence of coexisting monolayer phases |24j . The 
more general treatment of boundary motion applicable to 
a domain of any shape has, however, remained an open 
problem of considerable interest. 

In this paper, we adapt boundary integral techniques 
from fluid mechanics |2f| to the study of thin layers rest- 
ing on a subfiuid. Our approach is applicable to arbitrary 
geometries, and hence offers several advantages over tech- 
niques based on eigenfunction expansions that are useful 
only in particular situations with a high degree of sym- 
metry. It allows the comparatively straightforward cal- 
culation of linearized growth rates about any number of 
stable shapes; it is also a starting point for the detailed 
investigation of the boundary dynamics in the nonlinear 
regime. Further, by separating those aspects of the prob- 
lem that depend on energetics from those determined by 
the hydrodynamic equations, we gain several physical in- 
sights. Finally, our formulation provides a new example 
of a dynamical law governing curve motion in the plane. 
It is thus of interest to the broader study of pattern for- 
mation. 

In Section II, we present a formulation of the bound- 
ary dynamics valid for arbitrary viscosities and subfiuid 
depths, and note some global properties of the resulting 
dynamical law. Section III considers in greater detail the 
more analytically tractable limiting cases of an infinitely 
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deep or very thin subfluid. In Section IV, we examine 



the linear stability of straight and circular boundaries. 
Drawing on these results, Section [v] proposes a simplified 
dynamics appropriate when dissipation occurs mostly in 
the subfluid. This limiting form has the intriguing fea- 
ture that it contains a local drag term at the interface 
not unlike that introduced in the Rouse model of poly- 
mer dynamics [ p6| . Section VI gives a brief, illustrative 
calculation in the nonlinear regime, and Section VII dis- 
cusses some limitations and possible extensions of our 
work. 



II. GENERAL FORMULATION 

A. The Model and a Boundary Integral Formulation 

We begin by ignoring the presence of domains and 
studying a homogeneous two-dimensional layer coupled 
to a water subphase. The geometry of the system is 
shown in Fig. |l|. Both fluids are taken to have an in- 
finite horizontal extent, while the subphase may have a 
finite depth d. The monolayer is assumed to be an in- 
compressible, Newtonian fluid with surface viscosity r\ 
(dimensions mass/time) filling the plane z — pTfl . It 
rests on a three-dimensional incompressible fluid of vis- 
cosity rf (dimensions mass/ (length x time)) that occupies 
the region —d < z < 0. All variables referring to the sub- 
fluid are primed. Since the Reynolds numbers involved in 
the slow relaxation of micron-scale domains seldom ex- 
ceed 10 -4 , we are justified in working in the overdamped 
limit. Neglecting all inertial terms in the Navier-Stokes 
equation, we then find that the system is governed by the 
two coupled Stokes equations 



and 
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along with the incompressibility conditions V ■ u = 
and V • u' = 0. The term F s gives the body force 
that the subfluid exerts on the monolayer. Note that the 
monolayer is treated as a truly two-dimensional fluid, so 
p is a surface pressure with dimensions of force/length. 
The two fluids are also linked by no-slip boundary con- 
ditions; these simply require that, at z — 0, u' ± = u and 
u' • z = 0, where the _L indicates the in-plane component 
of a three-dimensional vector. In addition, we demand 
that all velocities vanish at infinity, and we impose no- 
slip boundary conditions on the bottom of the subfluid 
trough. Our model is thus essentially the same as that 
first introduced by Saffman to model flow in fluid mem- 
branes [^8| ; its predictions have been shown to agree well 
with experiment for monolayers flowing through a chan- 
nel JitJ . The form of the force F s is determined by the 
subfluid's Newtonian stress tensor evaluated at z = 0. 



The no-slip boundary conditions on the velocity field u' 
imply that all terms containing u' • z or its derivatives 
vanish, leaving the simple expression 



dz 



(3) 



The presence of a domain boundary will modify the 
equations for flow in a homogenous layer. We describe 
the interface as a closed curve C in the x—y plane (Fig. |l|) 
and ask that the fluid velocity be continuous across this 
boundary. In order to make the problem more tractable, 
we also assume that the two monolayer phases separated 
by the interface possess the same viscosity; this assump- 
tion will be discussed further later in this section. With 
the curve is associated a parametrization r(a) and an en- 
ergy £ [r] that is a functional of r(a). In the simplest case 
in which there exists only a line tension 7, for example, 
8 [r] = 7 Jds, where ds — y/gda, and yfg = \dr/da\ is the 
metric factor. 

To understand the effect of an interfacial force, we 
next introduce the Green's function that gives the 
response of a monolayer coupled to a subfluid to a point 
force exerted on the monolayer. A force g acting at the 
origin will induce a velocity field u s in the monolayer 
that is related to the Green's function by 
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(4) 



where summation over repeated indices is implied. One 
can similarly introduce Green's functions for the pressure 
and the body force exerted by the subfluid, defined by the 
relations 

P S M = i^iMft and = ■ (5) 

Together, the three will satisfy the equations 

V 2 G y - - diPj + /| = -47r%£(r) (6a) 



diG i:i = 



(6b) 



The first of these equations is the analog for Green's 
functions of equation (|l]) governing flow in the mono- 
layer, while the second reproduces an incompressibility 
constraint. One could of course also write down the equa- 
tions corresponding to the subfluid Stokes equation (^) 
and to the expression (|J) for F s , but they are not neces- 
sary for the further development of the present paper. 

For our purposes, the essential feature of an interface 
is that it exerts a force on the surrounding fluid. To 
find the velocity field in the presence of an interface, one 
must thus sum the contributions from the forces it exerts 
at each point. In other words, 



— J dsGaiT-Tia^AMTis)) , 
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where we have chosen the sign convention that Af is the 
net force per unit length that the fluid exerts on the in- 
terface. With equation (0), we have effectively reduced 
the problem of solving the two coupled partial differen- 
tial equations ([!]) and (^) to that of finding the correct 
Green's function; although equation (0) makes no ex- 
plicit reference to the subfluid, it is included implicitly 
because its presence determines the form of Gy. Vir- 
tual work arguments give an expression for Af in terms 
of £[r]; when a small element of the surface is moved 
slightly, the stresses in the surrounding molecular layer 
must do work to change £. The balance is expressed by 
the relation 
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(8) 



which is simply a generalization of the well-known Young- 
Laplace formula for the force exerted by a tense interface. 
The preceding two equations are the basis for all of our 
subsequent treatment of boundary motion. With them, 
we can calculate u anywhere in the plane; in particular, 
the values of u on C give the interfacial velocity. Observe 
that the boundary velocity is always determined nonlo- 
cally, and so depends on the shape of the entire domain. 

Equation ([?]) for the fluid velocity can also be derived 
by a more formal route j25[ : In the spirit of textbook so- 
lutions of Laplace's equation [ p9| , one begins by proving 
a reciprocal identity relating two arbitrary flows u and 
v. Choosing v to be proportional to G^, one finds that 
the limiting value of u as it approaches the boundary C 
of an arbitrary region is 
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ds Ui (r(s))T ijk (r-r(s))n k (r(s)) , (9) 



where fa is the force that the fluid inside C exerts on 
C, rik is a component of the unit normal vector, and 
Tijk = !•>,/, /', + dkGij + diGkj is the Green's function 
for the stress tensor. In the case in which the viscosi- 
ties inside and outside C are equal, one can readily com- 
bine the expressions for the limits from the inside and 
the outside to recover equation (Q). When the two vis- 
cosities are different, this approach does not generally 
work: Gy and T^k can depend on 77, so terms containing 
the Green's functions for the inside and outside regions 
cannot necessarily be combined and cancelled. Hence, 
except in certain limits, our theory cannot immediately 
be extended to include viscosity contrast. Because of the 
difficulty of measuring the viscosities of individual phases 
in the region of coexistence, the importance of such con- 
trast is usually not known. In most systems of interest, 
dissipation in the monolayer is negligible compared to 
dissipation in the bulk |lj,^|, suggesting that viscos- 
ity differences may not introduce too strong an effect. 
Likewise, in the opposite limit in which the subfluid is 



completely ignored, Gy does not depend on r\. The vis- 
cosity then enters the problem only through the factors 
of 1 jr\ in front of the integral in equation (||) , and a little 
algebra makes it possible to deal with variations in 77. It 
has been our experience that, in this case, the correction 
term due to viscosity contrast usually vanishes to linear 
order. Nonetheless, a more complete theory would allow 
for the possibility of viscosity differences between phases. 

Before turning to the calculation of Gy, we lastly re- 
mark that our curve evolution dynamic produces a gra- 
dient flow in configuration space. That is, the energy £ 
associated with an interface will always decrease mono- 
tonically in time, with d£/dt precisely given by the sum 
of the rates of viscous dissipation in the monolayer and in 
the subfluid. To prove this, we begin by observing that 
Af = Act ■ n, where Act is the difference in stress tensors 
across the interface. The time derivative of the energy 
is then d£/dt = Jds (6£/Sr) ■ u = J dsn ■ Act • u. The 
divergence theorem allows one to transform this integral 
into an integral over the plane; by using the dynamical 
equation ([!]), one can then make terms in F s appear that 
can by similar arguments be written in terms of integrals 
over the subfluid volume. Ultimately, one finds that 



^ = - 2?? /d 2 r(e l ,) 2 - -■>,,' 



d 3 r' U t 



(10) 



where = (diUj + djUi) /2 is proportional to the vis- 
cous part of the stress tensor in the monolayer, and a 
similar definition holds for e' in the subfluid. Summa- 
tion over repeated indices is implied. The two integrals 
give, respectively, the rates of viscous dissipation in the 
monolayer and in the subfluid. The result is hardly sur- 
prising, but it has potentially important consqences for 
the pattern forming properties of the model: The system 
is constrained always to move "downhill" in the space of 
shapes, so it will tend to get caught in metastable min- 
ima, and many shapes will be inaccessible to it. 



B. Calculation of the Green's Function 

To find Gij , we begin by finding the velocity field in the 
subfluid induced by an arbitrary flow in the monolayer; 
we will then proceed to calculate F s for this velocity field 
and finally to solve for the Green's function. The first 
task is greatly simplified by an observation of Stone and 
McConnell |13|, who showed that when the monolayer 
velocity field is incompressible the subfluid pressure p' is 
constant. Though their proof only holds for an infinite 
subfluid, it can readily be extended to the case where 
the depth is finite; the result can also be verified inde- 
pendently starting from the expression for the subfluid 
velocity as an integral over the plane z = p5| . With 
a constant pressure, Vp' = 0, and each component of u' 
becomes harmonic. Solving for u' is thus reduced to an 
exercise in electrostatics. Because the system is invariant 



3 



with respect to translations in the monolayer plane, the 
drag force must take the form 



FS M = |^ jd 2 vK{v -v)xi{v) , 



(11) 



where K is the derivative with respect to z of the ap- 
propriate Green's function for Laplace's equation. Both 
K and u are functions defined on the plane. When the 
depth d — > oo, K is readily obtained by the method of im- 
ages; the extension to the case of a finite depth is treated 
in Appendix |a|. It turns out to be convenient to pro- 
ceed in Fourier space. Denoting the Fourier transform 
(in two dimensions) of a function by a hat and adopting 
the convention that /(r) = Jd 2 qf(q) exp(iq • r), we find 
that 



K(q) 



■ coth(qd) 



whence by the convolution theorem 

F S (q) = -r)'qcoth(qd)u(q) . 



(12) 



(13) 



The obvious next step is to take the Fourier transform 
of the Green's function equation ([3a|). Making use of the 
fact that V • u = 0, one can obtain an expression for Pj 
and thus show that 



GiM) = 



irq 2 [q 2 + aqcoth(qd)] 



(14) 



The parameter a is the ratio of the viscosities of the sub- 
phase and the monolayer, 
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(15) 



and has dimensions of inverse length. This parame- 
ter plays a fundamental role in all subsequent analyses. 
When multiplied by an appropriate length scale, it will be 
the governing dimensionless parameter of the problem. 

The form of the expression for G y suggests that we 
define the differential operator 
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Dij = -- (%V 



(16) 



so that we can write G y - in terms of a single scalar func- 
tion as 



Ga = D l3 H(v) 



with 



#(q) 



q 2 [q 2 + aqcoth(qd)} 



(17) 



(18) 



In principle, we could now invert this transform and cal- 
culate Gij. In practice, the result would be so cumber- 
some as to be useless. Instead, in the next section we will 
study the behavior of H in several different limits. 



III. LIMITING CASES 

In this section, we will consider the behavior of the 
model in the limits of large and small a and d. Since 
both of these variables have dimensions, we must com- 
pare them with some other quantity to have a meaningful 
notion of "large" and "small." The only candidate that 
presents itself in the present formulation is the wavevec- 
tor q. If the system under consideration has a single 
length scale L, then the most important contributions 
in Fourier space are likely to come at q ~ l/L, and the 
relevant dimensionless parameters are d/L and aL. This 
would be the case, for example, for a circular domain 
of radius R whose boundary was smooth on all smaller 
scales; we would then have L ~ R. We caution, however, 
that q ~ l/L is at best a rough estimate and that q must 
actually be allowed to range from zero to infinity. Hence, 
great care must be exercised in taking limits when there is 
more than one candidate for L. This is particularly true 
when an approximation is valid only for large enough L, 
for there is always the danger that the domain boundary 
will finger or develop roughness at smaller length scales. 
For example, the important length scale in the case of 
a circle subject to n— fold harmonic perturbations is not 
its radius R but the wavelength 2tt R/n; whether or not 
a given approximation is valid thus depends on the mode 
one is considering. 



A. Infinite Subfluid 

We begin by considering the case of a very deep sub- 
fluid, qd — > oo. This is usually the experimentally rele- 
vant limit, for typical troughs have depths on the order of 
millimeters, while the monolayer domains observed tend 
to be on the scale of tens, or at most hundreds of mi- 
crons In this limit, coth(qd) — > 1, and H takes the 
simplified form 
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q 3 (q + a) 



(19) 



Note that in the limit of large q/a, H behaves as q~ 4 , 
while in the limit of small q/a it behaves as q~ 3 ; these 
should determine the behavior of Gy at small and large 
r, respectively. If one lets q/a — > oo, or equivalently sets 
a = 0, one recovers the case of a purely two-dimensional 
layer without any coupling to a bulk fluid. Although the 
integrals required to take the inverse transform of q~ 4 
diverge in the usual sense, they can be dealt with by the 
theory of generalized functions [^2|,^3). Essentially, all 
that is required is the introduction of a convergence factor 
like that commonly used to treat quantum-mechanical 
scattering from a Coulomb potential. One then finds 
that 



H(t) 



] ln(r) 



(20) 
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and 



G ij (r) = -8 ij hx(r) + -y- 



(21) 



where is a component of r. This is precisely the 
two-dimensional "Stokeslet" of fluid mechanics p5| , p4| . 
Although one should generally not take the logarithm 
of a quantity with dimensions, in the above equations 
this transgression turns out to be without consequences. 
With a = 0, the system has no intrinsic length scale, and 
we can chose to divide r by whatever length we please; 
the only change will be in an unimportant additive con- 
stant corresponding to a Galilean transformation. 

It is tempting to try to treat the limit in which dis- 
sipation in the subfluid dominates in the same manner 
just used for the case in which it is negligible. Unfor- 
tunately, if one blithely takes the inverse transform of 
l/(aq 3 ), obtaining — 27rr/a, and proceeds to calculate a 
Green's function, clearly unphysical results emerge. In 
the simplest terms, one can argue that this occurs be- 
cause the limit of large q/a will break down at large dis- 
tances, a feature that is acceptable for finite size domains, 
while the limit of small q/a becomes invalid precisely at 
the small distances that one must always consider. More 
specifically, the limit a — > oo corresponds to neglecting 
the term ?yV 2 u in the dynamical equation (|l|). One is 
then left with a lower order equation for which one is 
allowed to impose fewer boundary conditions, so the in- 
tegral formulation of equation (|7|) can no longer hold. 
We will argue in Section [v| that these difficulties can be 
circumvented with the judicious use of cut-offs. 

To proceed further within the present framework, how- 
ever, one must deal with the form for H valid for ar- 
bitrary a. A fairly involved expression for the inverse 
transform is obtained in Appendix [b|; it can be expanded 
about ar = to give 



H(t) = 



4C- 3 + 41n(ar/2) , n 
5 ( ar ) 



2 3 1 - C - Mar/2) 4 
■gW + 32 M + 



(22) 



where C ~ 0.577 is Euler's constant. Since Dijr 2 = 2Sij, 
the term proportional to r 2 will add a constant velocity 
to Gij . Unlike in the case a — 0, such a constant now has 
a physical meaning; the presence of the subfluid destroys 
Galilean invariance. 



B. Thin subfluid 

We next turn to the case of a very thin subfluid layer, 
qd — ► or d/L — ► 0. Although no experiments have yet 
been conducted in this regime, it seems plausible that it 
might be experimentally accessible. In numerical stud- 
ies of monolayer flow in canals, for example, Stone has 



observed that the effects of finite depth become impor- 
tant when d ~ L 18 1. For the largest experimentally 
accessible monolayer domains, a trough with a depth 
d ~ 100/im would then be required; this seems mechan- 
ically conceivable, although reflection from the bottom 
of the trough might make visualization with some mi- 
croscopy techniques difficult |l2). Beyond the fact that 
its theoretical treatment is less involved, this limit has the 
potential advantage that it would give the experimenter, 
in the depth d, an additional parameter that could be 
controlled with a fair degree of precision. For example, 
Klingcr and McConnell have reported the ability to set 
d to within 1/iin pl| . 

To lowest order in qd, H takes the form 



H 



q 2 (q 2 + A 2 ) 



A 2 = - = H- 
d i]d 



(23) 



The parameter A plays the same role as a in the infinite 
subfluid problem. The inverse transform of H can be 
taken without great difficulty and a Green's function ob- 
tained. It turns out to be more instructive, however, to 
back up several steps and to consider the function K in- 
troduced in the previous section. Expanding the Fourier 
transform of K, inverting, and taking the convolution 
with u, we obtain 



ird 



1 



3-M) 2 



K(r) 
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~d 



d 



6{r) + -d 2 V 2 S(r) 



1 



-d 2 V 2 u 



(24a) 



(24b) 



(24c) 



The first term in the expansion of F s has been derived 
by Stone via more heuristic arguments fl4| , |l8|| ; it is also 
the expression for F s that one obtains by treating the 
subfluid in the lubrication approximation. In qualitative 
terms, the series tells us that as d increases the force 
exerted by the subfluid develops an increasingly nonlo- 
cal character. For very small d, the no-slip boundary 
conditions on the bottom of the trough completely dom- 
inate behavior and prevent the effects of motion in the 
monolayer from propagating through the subfluid. As d 
increases, however, different parts of the adsorbed layer 
become more and more able to communicate with each 
other through the subfluid. In the opposite limit of in- 
finite depth, we thus expect that F s will depend on the 
velocity field throughout the monolayer. 

We confine ourselves for the moment to considering 
only the leading term in the above expansions. It is of 
course also possible to look at higher order approxima- 
tions, but they lack the internal consistency of the lowest 
order dynamics. In particular, they will not always yield 
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a law of motion that is a gradient flow, fn the first ap- 
proximation, a comparatively simple equation holds: 



?yV 2 u- Vp- — u = 
a 



(25) 



This equation has the same form as the Laplace trans- 
form with respect to time of the linearized Navier-Stokes 
equation pdu/dt = i]V 2 u — Vp and has been studied in 
this context [p5[ ; it has also been used to model flow 
in porous media |35[ |. If we ignore dissipation in the 
monolayer entirely compared with dissipation in the sub- 
fluid, the equation reduces to Darcy's law u cx — Vp, 
which describes quasi-two-dimensional flow in Hclc-Shaw 
cells [||J23). The Green's function for Eq. (§1 has previ- 
ously been calculated |25|] and can readily be obtained by 
the taking the inverse transform of equation ((2^). One 
finds that 



H(r) 



2n 



In (Xr) + K (Xr)} 



(26) 



and 



G y (r) = -26, 



+2- 



(Ar) 



— -Ao (Xr) 



(Ar-r 



K (Xr) 



Xr 

2Ki (Xr) 
Xr 



(27) 



where Ki is a modified Bessel function. Note that in the 
limit Xr — > 0, H approaches the expression obtained in 
the complete absence of a subfluid. Close enough to a 
singularity, the presence of the subfluid will always be 
negligible because the higher order derivative V 2 u will 
always dominate F s . 



IV. LINEAR STABILITY 

While the present formulation of the monolayer hy- 
drodynamics permits the investigation of boundary mo- 
tion far from any regular geometric shape, where a fully 
nonlinear analysis is necessary, we focus here on the cal- 
culation of linearized growth rates about several simple 
geometries. To develop physical understanding, we will 
begin with an interface in the shape of a straight line; we 
will then turn to the more experimentally relevant case 
of a circular boundary. 



A. Stability about a Line 

Consider an unperturbed interface that rests on the 
line y = and subject it to a sinusoidal perturbation 
y(x) — y(k)e ikx e y , where e y is a unit vector in the y- 
direction; we use the wavevector k to distinguish a one- 
dimensional Fourier transform from a two-dimensional 
transform with wavevector q. This gives rise to a force 



Af = f(k)e tkx e y which in turn causes a velocity at the 
interface u(k)e lkx e y , with u = dy/dt. For a boundary 
without internal structure, Af must always be normally 
directed, so there is no possibility that it will have an e x 
component. The relationship of / to y will depend on 
£[v], and its particular form is not of immediate interest 
in the present discussion. Quite generally, though, we 
expect that u oc / cx y. Since / is already first order 
small, once it has been calculated we may consider that 
the boundary takes its unperturbed shape. In those cases 
in which an analytic expression for Gij is known, we may 
then simply calculate the velocity component as 



u(k) 



f(k) 



dxG 2 2(xe x )e 



ikx 



(28) 
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We find that in the limit of negligible subfluid dissipation 
(o-fO), 



u(k) 



f(k) 
4ri\k\ 



while in the case of a thin subfluid, 



u(k) 



m 

2r/X 2 



W 2 



x 2 



(29) 



(30) 



Once again the effects of the subfluid are unimportant at 
small enough length scales: As \k\/X — > 00, the expres- 
sion approaches that valid in the absence of a subfluid. 
Similarly, as \k\/X — > 0, we recover the form that has pre- 
viously been calculated starting from Darcy's law ||22| , [23| . 

When a compact direct space form for Gij is not 
known, one may still find the growth rates by remaining 
in Fourier space. The method is presented in Appendix |c] 
for stability about a circle; the results of a similar cal- 
culation for a line are plotted in figure ^| in terms of a 
reduced growth rate a(k) = u(k)r/' / f (k) . At present, 
we simply state that, for an infinite subfluid in the limit 
\k\/a < 1, 



u(k) = 



irr/' 



(31) 



This is the same k dependence one would expect if dis- 
sipation occurred only at the boundary of the domain, 
instead of in the bulk fluid developments of Section 0. 

The k dependence of these growth rates can be under- 
stood on the basis of relatively simple arguments. Sup- 
pose that -ff(q) ~ q n ] n is determined by the number 
of derivatives of u and of p in the dynamical equation. 
Then, we expect that Gij ~ q n+2 . To find the linearized 
growth rate about a line, one must first take the inverse 
transform of Gij in two dimensions, then, in a rough 
sense, take a one-dimensional Fourier transform of the 
resulting function. On purely dimensional grounds, this 
will introduce an additional factor of q. Indentifying q 
with the wavevector k of the perturbation, we then ex- 
pect that u(k) ~ k n+3 . This is indeed the case: When 
a = 0, H ~ and u(k) ~ k^ 1 , and similarly for the 
other limits. 



G 



B. Stability about a Circle 

As in the previous section, we begin by considering 
a slightly perturbed domain parametrized as r(9) = 
R(l + e n e m6 )e r , where e r is a radially-directed unit vec- 
tor. We expect a force f n e m8 e r and a normal velocity 
component u n e e r ; u n — Rde n /dt. The force must 
again be normally-directed for a structureless interface. 
Incompressibility requires that the fluid velocity have a 
tangential component at the interface, unlike in the case 
of a line; this component does not, however, affect the 
evolution of the boundary's shape, so we will ignore it. 
To lowest order, we may still consider that the force acts 
at the unperturbed circle: Although there can now exist 
a zeroth order force, it must be independent of 9 and so 
will not cause any fluid motion, even when acting at the 
perturbed interface. 

Whereas for stability about a line we used a variety of 
direct space forms for G^-, here we shall instead derive 
one expression valid for arbitrary a and d. We begin by 
observing that, using polar coordinates and integrating 
over the polar angle, we may formally write 



H{r) = 2tt / 
Jo 



dq 



M ar ) 



q 2 [q + acoth(qd)] 



(32) 



Starting from this representation of H, straightforward 
but lengthy manipulations (described in Appendix |^) 
lead to an expression for the growth rate u n : 



n 2 Rf n 
V 



dw 



J 2 nH 



w 2 [w + aRcoth(wd/ R)] 



(33) 



of 



This is the central result of this section and the analo 
the growth rates given by equations (|29|), (|30|), and ( pi] ) 
for perturbations about a line. If we let the subfluid be 
infinite and set coth( wd/R) — 1, we recover the expres- 
sion previously obtained by Stone and McConnell [fL3| . 
When a — or aR — > oo, the integral can be evaluated 
in closed form. One finds that 



(H>2) 



4n n 2 — 1 

when dissipation occurs only in the monolayer, and 

4 n 2 .. 

fn (M > 2) , 



irrj' An 2 — 1 



(34) 



(35) 



when the subfluid dominates. The leading corrections to 
these expressions for finite aR can also easily be com- 
puted. In both cases, the first correction tends to de- 
crease the magnitude of the growth rate. This is not sur- 
prising: To zeroth order, we entirely ignored dissipation 
in the subfluid, in the one case, and in the monolayer, in 
the other. The next term, by accounting for these addi- 
tional sources of energy loss, increases the total amount 
of damping and so slows down the dynamics. 

One may verify that in the limit n — » oo, R — > oo, with 
k = n/R fixed, the results (|3) and (|35|) for a circle tend 



to the growth rates about a straight line ( |29| ) and (|31 
For n large enough, the curvature of the circular bound- 
ary is unimportant on the scales over which the forces and 
velocities vary, and the boundary acts essentially like a 
line. Hence, the length scale with which one must com- 
pare a is not the circle's radius R but the wavelength of 
the perturbation, which is proportional to 1/k = R/n. 
The expression (|3J) is the appropriate approximation for 
small aR/n, while ( |35| ) is more accurate for large aR/n. 
This dependence on n of the dominant source of dissi- 
pation can be seen in figure |^, where the exact reduced 
growth rate a n = u n rf / f n is compared with the two lim- 
iting forms. For fixed aR 3> 1, the growth rate is roughly 
independent of n for small n but decays like 1/n for large 
n. The crossover occurs where the curves given by equa- 
tions ( |34| ) and ( |35| ) intersect, at aR/n ~ 1. Except for 
a multiplicative factor, the growth rate about a line is a 
function only of a/k. This is not precisely the case for a 
circle, but, even for n small, many important quantities 
depend essentially only on aR/n. For example, figure |] 
plots versus aR/n the fractional difference between the 
exact growth rate given by equation (|3^) and the ap- 
proximation of equation ( pq ) for a number of different 
values of aR and n. To a good approximation, all of 
the points fall on the same curve. Actually, this collapse 
occurs over a much wider range of aR/n than shown in 
the figure, even when the fractional error is greater than 
0(1). Figure | also gives an idea of the error involved 
in using the approximate expression for the growth rate. 
We see that for aR as high as 100, the approximation is 
accurate to within 5% for the first few modes, with error 
increasing linearly as nj 1 aR for larger n. This asymptotic 
dependence of the error may also be derived directly by 
expanding the integrand of equation S3 ||36f . 



C. Example: Dipolar Forces 

In previous sections, we made no assumptions about 
the nature of the energy functional £ [r] associated with 
the boundary. Here we undertake a sample calculation 
for the functional of greatest experimental interest, that 
of a dipolar domain. In this case, we associate with the 
boundary not only a line tension 7, but also the electro- 
static energy stored in the electric field created by the 
dipoles. In reality, this electric field exerts a force not at 
the boundary but on the bulk dipolar fluid. Under the as- 
sumption that the domain has a constant dipole density, 
however, the electrostatic body force per unit volume can 
be written as — V</>, where <p is an appropriate poten- 
tial energy. After the introduction of a modified pressure 
Pm = P+4 1 PH i the equations describing the bulk flow are 
thus unchanged, and the electrostatic interactions only 
enter through their effect on the boundary conditions. 
One can show that this effect is correctly incorporated 
into our formalism if one simply views the electrostatic 
energy of the dipolar domain as a functional of its bound- 
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ary's parametrization. Several equivalent forms exist for 
the energy Ed of an arbitrarily shaped domain with (con- 
stant) dipole density /i per unit area . Of 

these, the most useful to us takes the form of the energy 
of interaction of two current loops: 

£ d [r] = -^- f dsi J ' ds 2 i(si) ■ t(s 2 ) 

x$(|r( Sl )-r(s 2 )|//i) . (36) 

Here C is the curve parametrized by r, t is the unit 
tangent to C, h is the thickness of the monolayer, and 
$(£) = sinli _1 (l/f) + £ - y/1 + t 2 . Adding to the elec- 
trostatic term the usual line tension energy jL, where L 
is the length of the curve, we arrive at an expression for 
£[r]. A fair amount of algebra then yields the force com- 
ponent /„. In the case of monolayers, the thickness is of 
molecular size, h ~ 10A, and the typical domain radius 
is R ~ lOjicm, so we are justified in taking the limit in 
which the aspect ratio p = 2R/h ^S> 1. It has then been 
shown that |L5| 



f = 1 
U ~ R 



in 2 - l) 



1 " 1 

-N B (1 - An 2 ) Y — ^ 



(37) 



Here the dipolar Bond number Nb = 2/i 2 /7 gives the 
relative importance of the electrostatic and line tension 
forces. Substituting this expression for f n into any of the 
growth rates calculated in the last section, one obtains 
a prediction for u n that can be compared directly with 
experiment. With values of the line tension on the order 
of 1 x 1CP 8 erg/cm p5|, a domain radius R ~ 50^m, and 
subfluid viscosity rf = 1 cp, we obtain interface velocities 
on the order of 1 /xm/sec for small Nb, and considerably 
less near the branching instability. 



small region around the point where we wish to know 
the boundary velocity, and it seems plausible that one 
might be able to neglect the variation of physical quanti- 
ties across this region. Then, the contribution transmit- 
ted through the monolayer itself to the velocity at a point 
r will be proportional to Af(r). In this picture, the ef- 
fective force at the interface is thus a consequence of the 
extremely small length scales over which dissipation in 
the monolayer is important. These allow us to take this 
dissipation to be essentially local compared with the dis- 
sipation in the subfluid which retains its very nonlocal 
character. 

To put these ideas into mathematical form, we begin 
by finding the limiting forms of Gij for small and large ar. 
These can be obtained by straightforward differentiation 
of the corresponding limits of H(r), derived in previous 
sections and in Appendix One finds that, for ar small, 



Gij behaves as 



- + C + In (or/2) 



7? ■ (38) 



while for ar large the appropriate expression is 



GfJr) 



(39) 



Note the important feature of Gfj that, unlike most of the 
Green's functions we have examined, it does not contain 
a term proportional to As a first approximation, 
we will suppose that there is a sharp transition between 
small and large ar behavior. That is, we will approximate 
the full Green's function as 



, , j G U V ) . r<v/a 

' ~ 1 G%(r) , r>v/a, 



(40) 



where v is a constant of order unity that will be deter- 
mined later. The expression for the velocity at a point r 
then becomes 



1 



V. DYNAMICAL LAW FOR MOTION 
DOMINATED BY THE SUBFLUID 

Two important features of the limit in which ah — > oo 
have already been emphasized: First, the monolayer vis- 
cosity represents a singular perturbation that must al- 
ways be taken into account near enough to a boundary 
or to a singularity. Second, linear stability results suggest 
that this limit can be partially understood in terms of an 
effective local dissipative force that opposes the boundary 
velocity at a given point on the interface. It is the pur- 
pose of this section to use these two observations to find a 
simplified boundary integral expression valid as aL — > oo. 
The physical ideas that will motivate the discussion are 
relatively straightforward: The Green's function Gij de- 
viates appreciably from its asymptotic large a form only 
when r 1/a. For a large enough, this describes a very 



Uj (r) . ...... 

47T77 J| r _ r ( s )|>,;/Q 



d S G£(r-r( S ))A/;(r( S )) 



1 



d3GfJr-r( S ))AfMs))- (41) 



4t"7 J\r-r(s)\<v/a 

This expression can be written in an alternative form 
in the usual case in which there are no tangential forces 
and one is only interested in the normal velocity. Then, 
the product G^ A/j is not singular as r(s) approaches r 
because the only nonvanishing term in the sum contains 
a factor of [(r(s) — r) ■ n] 2 in the numerator. Choosing 
the coordinate s = when r(s) = r and expanding for 
small s, one finds that r(s) — r ~ t — (l/2)«;ns 2 + • ■ • . 
Hence, [(r(s) — r) • n] 2 oc s 4 , cancelling the singularity in 
the denominator of G^- and making it possible to extend 
the integral of G^ to zero. For notational convenience, 
we pick axes such that n and t correspond to e^, and e y , 
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respectively, when s = 0. Then, the normal component 
of u at r is just u±, and we can write 

«.(•) = -j^jT«i<GS(r-r(«))A/,(r(«)) 



I 1 "/ 

47T7/ J| r _ r ( s )| <w / a 



Gfr(r-r( S )) 



-Gfi(r-r( S ))jA/ 4 (r( S )) . (42) 

Because Gf^ contains a factor of 1/a, the first inte- 
gral in the above expression will be proportional to 1/a. 
Thus, to leading order in large a, only terms of order 
1/a in the second integral are of interest. At this level of 
approximation, only constant and logarithmic terms in 
the integrand need be retained — that is, forces and ve- 
locities are taken to vary very slowly across the region of 
integration. Then, the second integral becomes 



v/a 



T- / dsAMr) 



1 

47T7/ J_ v/a 

A/i(i> 



ln(a|s|/2) 



hx{v/2) 



(43) 



where it should be emphasized that A/i (r) = Af (r) • e x 
is the force component at the point where we wish to cal- 
culate the velocity and so does not depend on s. Hence, 
to leading order in large a, the second integral yields a 
local force. We must still, however, specify the numerical 
value of the proportionality constant v[\n(v/2) + C— 1/4] 
in equation (fl3|), or, equivalently, of v itself. This can 
be done by demanding that the approximations of this 
section lead to the same linear stability results as the 
limits as a — > oo of the expressions valid for arbitrary 
a already obtained. Thus we ask that v be such that 
u{k) = —f(k)/%r)' (Eqn. ||l]) for stability about a line, 
and similarly for other geometries. This requirement de- 
termines the unique value v ~ 2.88. Likewise, with this 
choice of v, one can show that the proposed dynamical 
law preserves the area of domains and that the interfa- 
cial velocity of an arbitrarily shaped domain will vanish 
when the force is constant and normally-directed, both 
very desirable features. Substituting for v, we find the 
approximation to equation (]?]) that is the central result 
of this section: 



ITT]' 4TTT] J c 



r( S ))A/ 4 (r( S )) . 

(44) 



The first term can be associated with dissipation in the 
monolayer and is entirely local, while the second corre- 
sponds to dissipation in the subfluid. Note that the factor 
of 1/a in Gfj , when multiplied by the 1/f] in front of the 
integral, yields so the velocity does not depend at 

all on the monolayer viscosity. 

It remains to clarify when the results of this section 
are expected to hold — that is, what it means for a to be 



large. By considering the corrections of order 1/a 2 to the 
expression (^), one can convince oneself that these are 
negligible only when 



a 3> n and a 3> 



A/j(r) 
A/i(r) ' 



(45) 



where k is the curvature, A/{ = d(Af(r(s)) -e x )/ds, and 
both quantities are evaluated at r. These expressions 
are not surprising: For smooth interfaces, the curvature 
will be of the order of 1/L, where L is the size of the 
domain; Af[/Afx gives the length scale over which the 
force varies. As we have already pointed out, both are 
important length scales in the system. 

Unlike the Green's function valid for arbitrary a, Gfj 
has a compact analytic expression in direct space. The 
formula (^) thus presents considerable advantages for 
analytic and especially for numerical calculations. Such 
a simplification is all the more welcome because almost 
all experiments to date have been performed on sys- 
tems where the subfluid viscocity dominates. Admit- 
tedly, these improvements come at the expense of an 
approximation that is not perfectly controlled, but that 
nonetheless appears physically reasonable. 



VI. EVOLUTION OF THE PERIMETER 

Finally, we provide a simple illustration of the use of 
the present formulation for calculations in the nonlinear 
regime. We consider a domain with only line tension en- 
ergy, £ — jL, and calculate the time derivative of its 
perimeter L. This quantity is of particular interest be- 
cause L can be observed directly in experiments. 

Parametrizing the curve C in the standard way as r(a), 
and recalling that 5L/Sr(a) = —dt/da, we can write 



dL f dt 

— = -7 / da— ■ u(a) 

dt J c da 



. ^ j dsi J ds2n{s\)ni(sx) 

XGij(Ti2)K(s 2 )nj(s 2 ) , 



(46) 



where 1*12 = r(si) — r ( s 2), k is the curvature, and s is the 
arclength parameter. Since £ cx L, it will always be true 
that dL/dt < 0; the proof is identical to that already 
given for an arbitrary energy functional. The Green's 
function Gij can be any of the several forms we have 
already presented. One could equally use the approxi- 
mate form ( fl4| ) of the previous section for the velocities, 
finding that 



dL 
~dt 



7 



7"7 Jc 



dsn 



— — — / dsx ds 2 n(s 1 )n i (s 1 ) 
4*riJc Jc 



xG£-(r 12 )K(s 2 K(s2) , 



(47) 
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when the monolayer viscosity is small. These expressions 
for dL/dt are the analogs of the area-conserving "curve- 
shortening equation" previously derived under the as- 
sumption that dissipation occurs only at the domain 
boundary pl[ |. This simpler equation has already been 
fruitfully compared with experiment J^,|l5| , and we hope 
that our extensions will likewise prove useful in analyzing 
data. 



VII. DISCUSSION 
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Equation (|7|) encapsulates the focus of this paper: It 
provides a general prescription for finding the velocity 
of the boundary of an arbitrarily shaped monolayer do- 
main. We have discussed the form this equation takes in 
several limiting cases and have shown that it can be used 
to make quantitative predictions that can be compared 
with experiment. In particular, section ^ suggests an 
approach to the limit in which dissipation in the mono- 
layer itself is small; since most experiments have been 
conducted in this regime jl2|,[50| , its dynamics are of par- 
ticular interest. Our approximation to equation (Q) has 
the interesting feature that it contains a local term whose 
effect is equivalent to that of a drag force acting directly 
on the interface. From a theoretical perspective, it is 
clear that our somewhat ad hoc treatment of this limit 
leaves a number of open questions. For example, the is- 
sue of whether the dynamical law (|44|) gives precisely a 
gradient flow demands further elucidation, as does the 
application of our arguments to the case of a thin sub- 
fluid, which seems to reduce to Darcy's law when the sub- 
fluid viscosity dominates. In general, a more mathemati- 
cally rigorous treatment is desirable. These problems are 
currently under active study Q. The logical next step 
would be to undertake a detailed numerical investiga- 
tion of relaxation and pattern formation in the strongly 
nonlinear regime; techniques that have been developed 
for similar problems j?5) should make such studies possi- 
ble. To provide a complete test of the theory, it likewise 
seems useful to conduct experiments on systems in which 
the monolayer viscosity is sufficiently large that it can 
be measured independently |j^,[l3). Experiments carried 
out on a shallow trough might also be of interest. Finally, 
there are several physical effects that our theory does not 
pretend to include. Foremost among these are, first, the 
possibility that one of the monolayer phases might have 
a significant compressibility, and, second, the presence of 
thermal fluctuations in domain shape. The inclusion of 
thermal fluctuations, in particular, might have a signifi- 
cant qualitative effect on the pattern-forming properties 
of the model. We also have made no attempt to treat 
the dynamics of the various liquid condensed phases that 
show hexatic or other long-range order. Clearly, such 
order greatly increases the complexity of the problem. 



APPENDIX A: DETAILS OF THE DERIVATION 
OF THE GREEN'S FUNCTION 

This appendix gives a few of the steps that were omit- 
ted in the main text in the derivations of equations |ll| 
and |l2]. Since the components of u' are harmonic, they 
can be represented in the usual manner as 



u K r o) = J d 2 Y Ul {v) 



dG(r,r ) 

dz 



(Al) 



where Q is a Green's function for Laplace's equation that 
vanishes on the plane z = 0, and the integral is taken 
over this plane |2E(| . The boundary conditions are that 
u'j_ = u at z = 0, so the two dimensional velocity compo- 
nents Ui appear in the integrand. The Green's function 
Q is obtained by an eigenfunction expansion in Cartesian 
coordinates. With the hat denoting a Fourier transform 
with respect to the variables x — xq and y — jjq, one finds 
that 



A, s 1 sinh(g2; > )sinh[g( 

C/(q, z, z ) = . 

irq smh(qd) 



>< + d)} 



(A2) 



where z> and z< are respectively the larger and smaller of 
z and zq. Since F s = —r]'du'±/dz, comparison of equa- 
tions |ll] and Al indicates that K = d 2 Q /dzdzo, with the 
derivatives evaluated at z = zq = 0. The expression [l^ 
for K follows immediately. 



APPENDIX B: INVERSION OF H 

In this appendix, we obtain an analytic expression for 
the inverse transform of H — l/[q 3 (q + a)]. Begin by 
expanding in partial fractions: 



aq° a^q z 



1 



1 



a 3 q a 3 (q + a) 



(Bl) 



The inverse transforms of the first three terms are 
known |33| . To invert the fourth, we use the general re- 
sult that ^-^(q)] = -(l/r 2 )J : '- 1 [V2 5 (q)], where T~ x 
denotes an inverse Fourier transform and the derivatives 
in the Laplacian are taken with respect to q. We can 
then write 
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H(r) = % 



1 f°° x - 1 



(B2) 



Here, we have expressed the result in a dimensionless 
form by setting £ = ar and have performed the in- 
tegral over the polar angle in the inverse transform of 
\/q[l/(q + a)]. The integral in the fourth term converges 
in the usual sense and so can be evaluated numerically. 
It can also be expressed analytically in terms of a rather 
lengthy expression involving generalized hypergeometric 
functions p F s . These have known Taylor expansions that 
converge everywhere in the complex plane p4| . The ex- 
pansions lead directly to the series (|2^) . It can also be 
shown that to leading order for large £, H(r) = — 27r£/a 2 , 
as one would expect from limiting form of the Fourier 
transform for small q. 



APPENDIX C: LINEAR STABILITY 
CALCULATION FOR A CIRCLE 

We take as our starting point the integral expres- 
sion ( f{2| ) for H(r). Differentiating under the integral 
sign, we find that 



2?r / dq 
to 1 VI 



D tj J (qr) 



-\- a coth(qd)] 
dq 

o q 2 [q + acoth(qd)] 

rvr. 



(Cl) 



In order to find u n , we calculate the velocity at 6 = 0. 
Parametrizing the curve as r(6) = R(cos9e x + sin9e y ), 
we may write the velocity in the usual manner as the inte- 
gral around the curve of the Green's function multiplied 
by the force: 



Rfn 



/■2% 

/ dOGij (R [(1 — cos 9) e x — sin 9e t 
Jo 



(C2) 



where e r j is a component of e r . Substituting the integral 



expression (Cl) for Gij, setting w = qR, and performing 



some algebraic manipulations, one finds that 

a0cos(2n0J 



irrjR J q 2 [q + a coth(gd)] J 



wJ\{2w sin< 
2s\x\<b 



w cos 2 (/>Jo(2wsin< 



(C3) 



The identity 7rJ„(z) 2 = d(j)J (2zsm(j)) cosine/)) @ 
and two integrations by parts enable one to evaluate the 
integral with respect to <f>. The growth rate given in 



section |y| is then obtained by changing variables of in- 
tegration from q to w. The same approach can also be 
applied to other geometries. 
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FIG. 1. Schematic of the system studied, in which an in- 
terface C exists within a monolayer of viscosity n resting on 
top of a subphase with viscosity r\ and depth d. 

FIG. 2. Comparison of the reduced growth rates cr(fc) = 
u(k)n' / f(k) about a line in the limits of a very thin subfluid 
and of an infinite subfluid. Both curves are for a system with 
a — 10/xm . The dotted line was calculated using equation 
( ^p| ) with d = 10/xm; the solid curve was obtained by the 
methods of Appendix B in the limit d — > oo. The two curves 
approach each other and drop off like 1/k as k — > oo, but have 
markedly different behavior for values of 1/k on the order of 
typical domain length scales. 

FIG. 3. Reduced growth rate o n = u n rj / '/„ versus the 
mode number n for a circular domain with aR — 25 resting 
on an infinitely deep subfluid. The solid line gives the exact 
value, calculated from equation |33| The downward sloping 
dotted line gives the expression valid when dissipation in the 
subfluid is negligible (equation |34|) , while the other dotted line 
plots the expression valid in the opposite limit (equation fjEj ). 
The exact growth rate is always less than either of the approx- 
imations, but approaches the limiting expressions as in 2 
and n — » oo. 

FIG. 4. Fractional error [er„(a — » oo) — a n ]/a n versus aR/n 
for a circular domain resting on an infinite subfluid. Here 
a n = u n rj' I f„ is the reduced growth rate calculated from the 
exact expression (133 ), while a n (a — > oo) is obtained from the 
approximation (Bq) valid when dissipation in the monolayer 
can be neglected. The solid line gives the best-fit power law, 
[a n (a — > oo) - <7„]/cr„ ~ !A/(aR/n). 
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